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Abstract 

Using a spin-rotation invariant six-slave boson representation, we study the square 
lattice Hubbard model with nearest-neighbor hopping t and next-nearest neighbor 
hopping t'. We discuss the influence of t' on the charge and magnetic properties. In 
the hole-doped domain, we find that a negative t' strongly favors itinerant ferromag- 
netism over any incommensurate phase, especially in the strong coupling regime. For 
positive t' magnetic fluctuations are suppressed. A tight connection between frus- 
trated charge dynamics and large magnetic fluctuations is pointed out. A clear 

tendency towards striped charge ordering is found in the regime of large positive t'. 
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I. INTRODUCTION 



The spin rotation invariant slave boson representation is applied to the t — t' — [/-model. 
This model is expected to be relevant to the physics of high-temperature superconductors, 
since it includes a reasonable description of their band structure. It is a good candidate for 
the description of itinerant ferromagnetism too. Both behaviors are expected to occur in 
different regions of the phase diagram. Indeed it can be thought of as consisting of three 
characteristic regions, depending on whether the magnetic fluctuations are strong or not, 
and whether they are ferromagnetic or anti-ferromagnetic. The size of these regions is tuned 
by t'. In the non- interacting limit the role of t' is to shift the van Hove singularity that lies 
in the middle of the band for t' = to the lower band edge for t' = —t/2 or to the upper 
band edge for if — t/2 . The extension of this physics to the weak coupling regime has been 
extensively studied by Lin and Hirsch J^j, and Benard et al. p|, and Lavagna and Stemman 
0. They found that, for large negative t', the physics is dominated by strong ferromagnetic 
fluctuations in the low density domain, and by strong antiferromagnetic fluctuations in the 
vicinity of half-filling. Quantum Monte Carlo (QMC) simulations have been performed 
too. In particular Veilleux et al. |4j] confirmed this behavior, and thus put it on a stronger 
basis. They also established that the static and uniform magnetic susceptibility goes over 
a maximum when the system is doped off half-filling. Recently Hlubina et al. |J studied 
the same model at densities corresponding to the van Hove singularity, and found that the 
system is an itinerant ferromagnet for large negative t', and an itinerant antiferromagnet 
for small negative t'. Unfortunately these techniques can only be applied in the weak to 
intermediate coupling regime, because of the minus sign problem for QMC simulations, and 
because the RPA is intrinsically a weak coupling approach. For strong coupling one usually 
resorts to variational methods ||[7|], (for a recent discussion see ||). In order to cover the 
entire parameter range it is tempting to apply the Kotliar and Ruckenstein slave boson 
approach . It not only proved to yield ground state energies very close to the exact 
ones, but very realistic values for the structure factors too JTIJ. Until now little attention 
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has been paid to the charge structure factor. The aim of this brief report is two-fold. 
First we calculate the g-dependent magnetic susceptibility in order to determine in which 
domain of the phase diagram the antiferromagnetic, ferromagnetic and incommensurate 
fluctuations dominate, both in the intermediate and strong coupling regime. Second we 
calculate the charge structure factor. We then show that strong magnetic fluctuations are 
systematically accompanied with a clear reduction of the charge structure factor, i.e. by a 
frustration of the charge dynamics. In this model this happens for negative t', in the hole 
doped region. For positive t', the charge structure factor is enhanced and the magnetic 
fluctuations suppressed. We obtained this result using the spin-rotation invariant (SRI) 



slave boson representation of the Hubbard model fll| , |T2|| , and the expression for the spin 
and charge dynamical susceptibilities we recently derived and applied to the Hubbard model 
|10| . We note that this representation has been recently revisited by Ziegler et al. |13 . Our 



expressions for the susceptibilities remain unchanged by their considerations. 

II. FORMALISM 

In this work we calculate the spin and charge dynamical susceptibilities of the two di- 
mensional t — t' — U model. The Hamiltonian reads: 

H = - J2 kj4aCja + n A n H (!) 

i,j,a i 

We consider the case where the hopping integral is tij = t for nearest neighbors, tij = t' 
for next-nearest neighbors and tij = otherwise. To this aim we apply the spin-rotation 
invariant (SRI) slave boson formulation of the Hubbard model [P1 , |l"2"f to the one-loop calcu- 
lation of the susceptibilities that we applied to the Hubbard model. In this framework the 
dynamical spin susceptibility is given by 

Xs\K,Uj) = =; =; =; , (2) 

1 + A kXo(k, w) + A lX i(k, u) + A 2 [xi(k, w) - Xo(K u)x2(k, w)] 

where 
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Xn(k, u) = - Yl (tp + tp+kTGoAPi iuJ n )G a(p + k,uj + iu n ) (n = 0, 1, 2) (3) 

p,iu> n ,o 

In the low frequency regime Eq. (0) has an RPA form, to which it reduces in the weak 
coupling limit. It nevertheless differs from it in two important respects. First the effective 
interaction Aj: does not grow indefinitely as U grows, but saturates at a fraction of the 
average kinetic energy in the strong coupling regime. Second it is fc-dependent. Therefore 
if an magnetic instability of the paramagnetic phase at a given density develops towards an 
incommensurate phase, characterized by a wave vector q, this wave vector will be different 
from the wave vector p at which Xo(p, 0) reaches its maximum. Thus at a given density, 
the wave-vector q characterizing the phase towards which the paramagnetic phase can be 
unstable to, depends on the interaction strength, in contrast to the ordinary RPA. It nu- 
merically turns out that the fc-dependence of is enhanced by increasing the interaction 
strength. A^ is typically largest for k = 0, and such is Xi{k)- The contribution involving 
A 2 is smallest, and has little influence on the magnetic properties. The numerous undefined 



symbols in Eqs. (0-0) can be gathered from Ref. |I0 |, except for a misprint there: the third 
line of Eq. (A8) should read: 

9^ = 2VW d + x+ te*\ (4) 



III. RESULTS 

We now proceed to the numerical results. We first calculate the density dependence of 
the static (but g-dependent) magnetic susceptibility. In order to magnify the effect of t', 
we perform the calculation for t' = — 0.47£. Had we chosen t' = —0.5t, then the van Hove 
singularity would lie right at the lower band edge. For U = At and (3 = 2 we display the 
density-dependence of Xs for several ^-vectors in Fig. [I]. For these parameters the paramag- 
netic phase does not show magnetic instability. At a particular doping the maximum of Xs 
(in its ^-dependence) tells us towards which phase an instability will develop. We checked 
numerically that this really happens at lower temperature. In the vicinity of half-filling Xs 
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is largest for the commensurate vector Q = (n,ir). In the low-density range Xs is maximal 
for q = 0. This range is very large and extends from S ~ 0.38 to S = 1, S being the hole dop- 
ing. In this domain the fluctuations are predominantly ferromagnetic, because the system is 
making use of the van Hove singularity to reduce its free energy. Between these two regimes 
there is a small window where the instability is towards an incommensurate phase with q 
along the diagonal of the Brillouin zone. We find that there is a value of the doping, that 
we denote So, beyond which Xs is largest for q = 0. So is seen to decrease with increasing 
interaction. We note that the doping So at which Xs is largest for q = decreases with 
increasing interaction. For U — we found it to be So = 0.42, while for U > 15t Sq goes to 
zero. We thus obtain that, for strong coupling, the paramagnetic phase is unstable towards 
ferromagnetism over the entire doping range. This dependence of So on U can be traced back 
to the g-dependence of the effective interaction entering Eq. (0), as discussed below Eq. 
(H). It turns out that the g-dependence is weak for weak coupling, and gets stronger with in- 
creasing U. This plays a crucial role in assessing towards which phase a magnetic instability 
may develop. We note that this effect is neglected in the usual RPA and in the two-particle 
self-consistent approach jTJ]. In those approaches all what matters is the g-dependence of 



the bare susceptibility Xo- We note that including t' changes dramatically the phase diagram 
as compared to the t' = case. In the latter case ferromagnetism may only show up for 
very strong coupling (U > 66i) and in a narrow doping region located around S ~ 15% |15| . 



The influence of i! on the doping dependence of the uniform susceptibility is displayed on 
Fig. for U = 4t and Fig. |3| for U = 20t. For moderate coupling decreasing t' changes the 
monotonic behavior of Xs{S) into a non-monotonic one, which is typical of high-T c materi- 
als. The height of the maximum increases with t', and its location is shifted towards higher 
doping. This behavior, as well as the location and the height of the maximum, agree with 
the QMC data of Veilleux et al. [|J. We thus conclude that the non-monotonic behavior 
of Xs for moderate coupling mostly results from band-structure effects. Experimentally the 
non-monotonic behavior of xs has been observed in La2_ x Sr x Cu04 and the maximum is 



reached for x ~ 0.25 ||16|| : in YE^CusOr-^ xs only increases with increasing hole doping, 



and one may assume that a maximum is reached for doping values that cannot be reached 



experimentally. According to Hybertsen et al. [[17]], t' = — 0.16t is relevant to La2_ x Sr a .Cu04, 
in which case our approach yields the location of the maximum of xs at 5 ~ 0.20. However 
the dependence of xs on 5 is too weak to reproduce the experimental data. In the strong 
coupling regime Xs has a maximum in its doping-dependence for t' — |T3]. Decreasing 



t' results into an enhancement of the maximum of Xs, and into a shift of it towards larger 
doping. For large t' its location coincides with the van Hove singularity. Accordingly this 
non-monotonic behavior of Xs for strong coupling results from a combination of interaction 
and band structure effects. We calculated Xs(o) f° r other values of q too. For t' = —0A7t 
it turned out that Xs(q = 0) is largest over the entire doping range. Thus raising up the 
interaction leads to a dramatic widening of the ferromagnetic domain. This is in agreement 
with the variational calculation of Pieri et al. ||, who investigated in more detail the low 
density route to ferromagnetism due to Muller-Hartmann J/]]. 

We now turn to the charge structure factor. The latter is calculated according to Eq.(13) 
and Eq.(17) of Ref. fllPf . We recall that the comparison of the slave boson charge structure 



factor to existing QMC data displays a quantitative agreement. Here we perform the cal- 
culation for finite t', for U = At, (5 — 8 and quarter-filling and display the result on Fig. 

As compared to the t' = result, decreasing t' (i.e. t'/t becoming increasingly nega- 
tive) substantially suppresses the charge structure factor, especially around (0,7r), but also 
around (n, it) for the largest t'. This suppression results from a frustration of the charge 
dynamics. We note that this suppression takes place in the parameter regime where the 
tendency towards magnetism is strongest. This leads us to propose that this is a general sit- 
uation: frustrated charge dynamics and magnetic instabilities are occurring simultaneously 
in strongly correlated systems. This is well known for the Hubbard model at half-filling 
where the charge dynamics is so strongly frustrated that the system becomes insulating and 
the physics is dominated by strong anti-ferromagnetic fluctuations. Here the effect is less 
dramatic since the system remains metallic, but clearly noticeable. In the opposite case of 
a positive t' (i.e. t'/t > 0) the charge structure factor gets essentially flat along the side 



of the Brillouin zone. It is particularly enhanced in the vicinity (0,7r). This can be partly 
understood using a local picture. Assuming that the charges can form either a checkerboard 
or a stripe pattern, one sees that the number of frustrated bonds is larger for the former 
pattern. This plays little role if t' is negative, since one hop along the diagonal costs en- 
ergy, but an increasingly important one when t' is increasingly positive. Clearly no such 
charge ordering occurs here, but tendencies towards such patterns emerge. The frustration 
effect is absent, and we indeed did not find any sign of a magnetic instability. The correla- 
tion frustrated charge dynamics-enhanced magnetic fluctuations can be traced back to the 
relationship between the local susceptibilities and the density 



which follows from the Pauli principle. Indeed if the local charge susceptibility is reduced, 
the magnetic one is enhanced, and vice versa. In Fig. [5] we display the density dependence 
of the charge structure factor, for U = 4t, (3 = 8 and t' = —0A7t. Under an increase of the 
density, S c first goes up, until 5 ~ 0.25 and then goes down while going closer to half-filling. 
In weak coupling one would expect S c to decrease upon doping, but the opposite behavior 
holds in a dense strongly correlated system. 

In summary we studied the charge and magnetic properties of the 2-dimensional t—t' — U 
model. We found that a negative t' has a strong influence on the phase diagram, a very 
large portion of it being dominated by strong ferromagnetic fluctuations. We also showed 
that frustrated charge dynamics and strong magnetic fluctuations occur simultaneously. We 
found a tendency toward striped phases only for large positive t'. 
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FIGURES 
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FIG. 1. Doping-dependence of the magnetic susceptibility for U = 4t, (3 = 2 and t' = — 0.47t 
for a set of wave- vectors lying along the three main directions of the Brillouin zone. Inset: Dop- 
ing-dependence of the magnetic susceptibility for U = for the anti-ferro- and ferromagnetic 
wave- vector. 
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FIG. 2. Doping-dependence of the static and uniform magnetic susceptibility for U = 4t, 
(3 = 3 and t' = -0.lt, t' = -0.3t and t' = -0.471 
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FIG. 3. Doping-dependence of the static and uniform magnetic susceptibility for U = 20t, 
(3 = 2 and t' = -0.lt, t' = -0.3* and t' = -0.47*. 
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FIG. 5. Charge structure factor of the t-t'-U model for U = At, (3 = 8 and t' = -0.47t 
several values of 5. 
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